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We study the decoherence of a quantum computer in an environment which is inherently correlated 
in time and space. We first derive the nonunitary time evolution of the computer and environment 
in the presence of a stabilizer error correction code, providing a general way to quantify decoherence 
for a quantum computer. The general theory is then applied to the spin-boson model. Our results 
demonstrate that effects of long-range correlations can be systematically reduced by small changes 
in the error correction codes. 
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Quantum computers bear the promise to solve certain 
problems exponentially faster than their classical coun- 
terparts 0. Although small computers have been suc- 
cessfully tested 0, the development of large computers 
has been hindered by decoherence. The most promising 
method to tame decoherence is quantum error correc- 
tion (QEC) 0000. In QEC, it is usually assumed that 
correlations in the environment cither are non-existent or 
decay exponentially in time and space. In contrast, re- 
cent work argues that correlated environments can lead 
to quadratically worse error- levels 0000- Though 
the assumption of uncorrelated noise is often reasonable, 
it is not fulfilled in several physical systems proposed for 
realizing quantum computers, notably solid state systems 
using superconductors or quantum dots [I(tI |. Hence, 
it is far from clear how much protection from decoherence 
QEC gives in these important cases . 

In this paper, we consider the long time dynamics of 
a quantum computer immersed in a correlated quantum 
environment and protected by QEC. First, we describe 
the parameters which quantify the level of protection 
from QEC and give explicit expressions for them. Sec- 
ond, we calculate these quantities in a concrete example: 
the spin-boson model [ll| . This model is directly applica- 
ble to solid state quantum computers 0E3 but formally 
outside the scope of QEC 0. 

Our work shows that some protection against long- 
range correlations can be built into QEC codes. The 
new element here is that the periodic measurements in 
the QEC method separate the environmental modes into 
high and low frequencies. This natural "new" scale can 
then be used to engineer quantum codes to better cope 
with the long-range correlations. 

To follow the long time behavior of the computer, we 
remove non-essential elements and assume: (1) Quantum 
gates are perfect and operate much more quickly than the 
characteristic response of the environment. (2) States 
of the computer can be prepared with no errors. (3) 
Thermal fluctuations are suppressed. Finally, for clarity 
in the spin-boson example, we consider ohmic coupling 
between the environment and the qubits. Extensions to 
sub-ohmic and super-ohmic coupling are straightforward. 



Decoherence, QEC, and correlations — The wave func- 
tions of the computer and the environment are unavoid- 
ably entangled during their time evolution. When a mea- 
surement is performed, this entanglement is translated 
into the probability of an error with respect to the ideal 
state of the computer. This general feature of open quan- 
tum systems is known as decoherence. In order to reduce 
its effects, QEC encodes the computer state in a subspace 
of a larger Hilbert space, TL, such that by measuring ob- 
servables in Tl the wave functions of the computer and 
the environment are disentangled to some order in e, the 
probability of a single error in a QEC cycle. The outcome 
of the measurement, called the syndrome, discriminates 
among the possible errors in the computer, allowing a 
correction to be made. Note that by performing a mea- 
surement on the computer, the environment's dynamics 
is restricted; thus, it is crucial to follow the non-unitary 
evolution of the entire system. 

The difficulty with correlated environments is that e 
is no longer time independent 0000: past errors 
or gates can substantially change e(t). Consider, for in- 
stance, that at time t\ an error occurs in the computer. 
Later, at time £2, this is realized and corrected. Although 
the error was detected, its precise time is unknown; there- 
fore, when calculating the probability of a future error, all 
possible events at t < £2 must be taken into account. If 
the environmental correlations decay algebraically, there 
is no characteristic scale at which to truncate this sum. 
Hence, it may be hard to estimate e(t). Sometimes it is 
sufficient to evaluate upper bounds on the error ; how- 
ever, these estimates may be poor for certain models. For 
instance, correlated noise was studied in 0, which de- 
rived an upper bound on the error strength in the Ohmic 
spin-boson model. This bound is linear in the bosonic ul- 
traviolet cut-off, A, which is the only scale in the model. 
One of our main points is that the dynamics imposed by 
the QEC code provides another scale, which can then be 
used to generate better codes and/or bounds. We specif- 
ically discuss the spin-boson case, but our results can be 
adapted to other situations. 

Time evolution of encoded qubits — The environment 
usually can be described by a Hamiltonian, Hq. Although 
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other couplings between the computer and the environ- 
ment are possible, we focus on the case of local vectorial 
coupling: V = ^ J2 X f ( x ) ' ( x ) > where a (x) are Pauli 
matrices for the qubits, A is the coupling strength, and 
/ (x) is a function of environment operators. 

QEC is essentially perturbative in V. It is therefore 
natural to define the interaction representation: V (t) — 
en Hot Ve~ ~R Hot . In this representation, the time evolution 
of the system during a QEC cycle that starts at time t = 
and ends at time t = A is 



U(A,0)=T t e 



(1) 



where T t is the time ordering operation. It is also simple 
to describe gates that are done faster than the environ- 
mental response. If a gate R T is executed at time r < A, 
at the end of the QEC cycle the time evolution is given 
by U(A,t) R t U{t, 0). 

At time t — A, the measurement selects a set of 
terms from the r.h.s. of Eq. (JTJ. Only this set is car- 
ried to the next QEC cycle. For stabilizer codes Q, it 
is straightforward to identify which operators must be 
kept. First, because the measurements are performed 
on individual logical qubits, it is sufficient to analyze 
the n physical qubits that define a logical one. All uni- 
tary operations in the Hilbert space of these n qubits 
can be written using the Pauli group G n . The sub- 
group E of G n of all possible errors in the computer 
is called the error set. Intuitively, E is given by the 
Pauli matrices appearing in each term of an expansion 
of the the r.h.s. of Eq. (JTJ. For instance, the element 
g = a x (xi) a z (x 2 ) € E is given by the time-ordered term 

A " C d< l Jo 1 dt *f X ( X L *0 f Z ( X 2>*2) O* (Xl) Z (X 2 ). 

The next step is to decompose E according to the pos- 
sible values of the syndrome. This is done by noticing 
that the measurement will not distinguish elements in- 
side the logical Hilbert space. Thus, the logical Pauli 
Group Gi, generated by the identity and the logical 
Pauli matrices, {l,X,Y,Z}, defines the subgroup E = 
G\ n E. The partition P of E given by all left cosets 
of Eq in E sorts the errors by their syndromes. Us- 
ing this fact, the r.h.s. terms of Eq. (JTJ can be re- 
ordered as U (A, 0) = J2 m u m (A, 0) with respect to 
P. Finally, a QEC code defines an appropriate recov- 
ery operation for each syndrome, R = {r m }. At the 
time of the measurement, t = A, only one {u m } is se- 
lected and the corresponding operation, r m € R, is per- 
formed. The overall time evolution for a QEC cycle is 
v m (A, 0) = r m (A) u m (A, 0); thus, interference between 
terms with different syndromes is eliminated. 

The analogous result for many logical qubits and a se- 
quence of N QEC cycles follows directly from the above: 



T 



v WN (NA,(N-l)A)...v Wl (A,0) 



(2) 



where w is the particular history of syndromes for all 
the qubits and v w . is simply the product of evolutions 



for the individual qubits. Each history comes with the 
associated probability 

P(T w ) = (poK^|T+ r T w 1^)1^,) (3) 

where tpo and ipo are the initial states of, respectively, the 
environment and the encoded qubits. Finally, QEC only 
partially disentangles the environment and the qubits. 
There is always some residual decoherence, which can be 
found from the reduced density matrix 

(^o|[^o|Tt|s)(r1T w |^o)] \<po) 
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with r,s elements of the logical subspace. 

Eqs. ©-(0J, our central formal result, can be used to 
assess the protection offered by a QEC code. They show 
that A is a natural scale in the field theory that describes 
the system's evolution. Although the consequences of 
this imposed scale are model dependent, its existence can 
be used to construct better QEC codes; to show this 
explicitly, we now turn to an example. 

Decoherence in the spin-boson model — The spin-boson 
model deals with generic two state systems coupled lin- 
early with an infinite set of harmonic oscillators. In its 
general form, the model includes a finite tunneling am- 
plitude and a bias between the two states. For qubits, 
these two features model "imperfections" and so are not 
considered here; they are not fundamental to the under- 
standing of decoherence in a correlated environment |l2l | . 
As we consider only "perfect" qubits, the computer expe- 
riences errors only due to dephasing. Finally, we choose 
to consider the case of linear coupling to an ohmic bath. 
We stress that these two choices do not restrict our re- 
sults, but allow for a convenient notation. 

The Hamiltonian of the model can be written as 



H = V ± dx[d x( j>{x)f+[Jl (x)f- 

^ J-oo 

(5) 

where <fi and II = d x 9 are canonical conjugate variables, 
cr^ act in the Hilbert space of the qubits, is the ve- 
locity of the bosonic excitations, and K = ks = 1. The 
bosonic modes have an ultraviolet cut-off, A, that defines 
the short-time scale of the field theory, t uv = (Avb) . 
Following our general assumption, we regard gates as per- 
fect and with operation time t g <C t uv . 

Between gates, the exact time evolution of a qubit in 
the interaction picture can be expressed as the prod- 
uct of two vertex operators of the free bosonic theory, 



U n (t, 0) 



T T \[8(n,t)-6{n,0)]c 



', and a coherent evolu- 



tion that is irrelevant to our discussion |12(. Hence, it 
is straightforward to express the reduced density matrix 
as a bosonic correlation function. In particular, for a 
single qubit in the initial state |-0o) = «|T) + P II), the 
off-diagonal element is 



(A) = a/3* (0| e V2^A[0(n,A)-0(n,O)] ^ = ^ 



(6) 
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FIG. 1: Steane's 3 qubit quantum error correction (QEC) 
code Q- The initial wave function, \i/j ) <8> (IT) + |l})/2 ® 
(|T)+|l>)/2, is encoded by two controlled-NOT (CNOT) gates, 
i?CNOT = af cr]+af (T^ , into an entangled state Vcncodc) = 
a |f) +0 |I) with |f) = (|fff) + |TU> + Iff I) + lilt)) /2 and 
|I) = (111!) + ITT) + |Ut) + ITU)) /2- After some time, the 
information is decoded by a second pair of CNOT gates. An 
error in is identified by measuring the the value of erf and 
erf (rectangle). The cycle of QEC ends with the correction of 
a possible phase-flip (arrow). 



where |0) is the bosonic vacuum, e = A 2 hi[l+(A Wb A) 2 ]/2 
is the probability of an error, and Ai^A ^> 1. 

QEC and the spin-boson model — The simplest QEC 
protocol is Steane's three qubit code (see Fig. QJ. It is 
designed to protect a logical qubit from a single phase- 
flip, that is from dephasing in lowest order in the coupling 
to the environment. Thus, we illustrate our discussion of 
QEC in correlated environments by applying this code to 
three qubits that dephase according to Eq. J5J. 

The error set of a single logical qubit is E = 
{I, a), o)al^ afalal}, where j, k = {1, 2, 3}. The log- 
ical Pauli matrix Z — cfa^i is contained in E. Thus, 
the subgroup Eq — {/, can be used to partition the 
error set into four equivalence classes, each with its re- 
spective recovery operation, P <-> R: 



{I, Z}^I , K, CTfcrf} <- <J» 



W2 



1} 



(7) 



The first class, {/, Z}, corresponds to a superposition 
of the state where no error occurred, (^encode), and that 
with a logical phase flip, Z IVencodc)- Similarly, the other 
three classes are superpositions of states with one and 
two phase-flips on physical qubits. When the recovery 
operations are done, the logical qubit is once again a 
superposition of IV'cncodc) and Z IV'cncodc)- Although the 
final state is seemingly the same, each possible evolution 
has a different bosonic content. 

The bosonic parts are written compactly by sep- 
arating the time evolution operator into "even" 
and "odd" orders: C/j(A,0) = rjj + iVjPj, 
where r, 3 = cos [y^A [9 (j,t) — 0(j,0)]] , v 3 = 
sin [yf \[6(j,t)-9(j,0)]]. Thus, the time evolu- 
tion by the end of a QEC cycle is 

w (A, 0) = ?7i7?2?73-f - iviV2V 3 Z, 

ui(A,0) = iv^q^l - rjiviVsZ, (8) 

Wa(A,0) = i-q^m 1 ~ v \ or 

f 3 (A, 0) = £771 7721^ - VxWNzZ. 

Given a particular history of syndromes, w, it is 



straightforward to find both its likelihood, V (T w ), and 
the corresponding contribution to the off-diagonal ele- 
ment of the reduced density matrix, p^j_ (T w ). For the 
moment, we consider the case of qubits separated by a 
distance much larger than i) D iVA, so that spatial corre- 
lations can be disregarded. Although a very particular 
example, this is illustrative of the general discussion. 

By the end of the first QEC cycle, there are two differ- 
ent cases to consider: T w = vq with probability V (vq) ~ 
1 — 3e/2 and T w = ^(1.2,3) with V (v n .2,3.)) — e/2. Each 
one causes an intrinsic dephasing for the logical qubit 

3p- e -t- P - 3e 

Pn (v ) = a(3* 



i a(3* (1 
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1 + 3e" 2e 

PTi (v ( i,a,3)) = ap*e- e ~ ap* (1 - e + ...) . (9) 

Note that by the end of a cycle with an "error", the 
logical qubit has suffered precisely the same decoherence 
as for an unprotected physical qubit, Eq. ©. A cycle 
with such decoherence has, however, become a "rare" 
event; this exemplifies the benefit provided by QEC. 

If there were no correlations (or in this case mem- 
ory) between QEC cycles, then the likelihood of an error 
would be history independent and given by @. In this 
uncorrelated limit, fault-tolerance can then be proved 

m 

Correlation between QEC cycles — We now address to 
what extent correlation between cycles changes the likeli- 
hood of errors. In the current example, the contribution 
to V (T w ) due to long-range correlation is elegantly eval- 
uated by the Operator Product Expansion (OPE) [l3j . 
The idea is to express the evolution of the environment 
during each cycle as a single operator that will capture 
the most relevant contributions to the correlation func- 
tions that need to be evaluated. 

In order to use the OPE, V (T w ) must be re-ordered 
such that operators with arguments closer in space and 
time are put together. Hence, in the spin boson case, all 
local operators are products of rjj and The OPE de- 
composes rjj and into high frequency, A -1 < uj < t>bA, 
and low frequency, u> < A -1 , parts. The high frequencies 
are integrated out and the most relevant components of 
the low frequencies are kept: 



lie" 



:cos 



27rA9 t 0(j,O)A]:}/2, (10) 



where : : 
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represents normal ordering, A + = rjj and A_ = 
In essence, this procedure captures the long range 
effects on the environment of the dynamics imposed by 
the QEC code by using operators which have the "new" 
ultraviolet cut-off A -1 . The leading order effects of cor- 
relations in V (T w ) are given by the effective operators 



1 - 3e/2- 
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V 3={ 1,2,3} 
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[d t eu,o)Y 



(ii) 
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Hence, the probability, P(T W ), of a history w is eas- 
ily evaluated using Wick's theorem and the fact that 
{:[a t e( J ,i)f::^( J ,0)f:>^l/(2,V). 

The simplest case is to calculate the probability of 
having errors in the QEC cycles starting at times t\ 
and t%. From Eqs. (fill) and (|12f> . it is straightforward 
that, in leading order in A, V (...Vj....Vj...) w (e/2) 2 + 
A 4 A 4 /[8(ti — t2) A ], where the first term is the uncorre- 
cted probability and the second is due to correlations be- 
tween errors in different cycles. Therefore, for N -C 1/A 2 
QEC cycles, the probability of having two errors of any 
kind is V 2 ~ [(e/2) 2 ^ + ^-N]. Thus, for an ohmic 
bath and finite operation time, correlations give a small 
correction to the usual distribution of probabilities. 

Reducing the effects of long-range correlations — It is 
well known that decoherence of a physical qubit can be 
systematically reduced by applying a series of NOT gates 
at a frequency higher than the bosonic cut-off 0|. Al- 
though possible, the restriction ui WbA can be experi- 
mentally very stringent. The analysis here implies that 
essentially the same effect can be obtained for a logical 
qubit with VbA replaced by A -1 . For instance, a simple 
logical NOT, X — erf erf erf, executed at half of each QEC 
cycle will change dramatically the effects of correlations 
in Steane's 3-qubit code. Following the same steps as 
above, we now find that although the local probability of 
an error increases to ~ 3e/2, the operator that captures 
the correlation between cycles is a second derivative of 
the bosonic field. Hence, Eq. (|12|l would now read 

«*={1A3} ~ 3e / 2 + (^A 2 A 4 /32) : [8 2 6] 2 :, (13) 

with Vq = 1 — Ylj={i 2 3} u j • The relevant two-point 
correlation function now decays faster, ~ 1/i 8 , implying 
a much reduced effect of correlations. 

In the general case of an environment with spectral 
function J(oj) ~ uj s , this simple change in the QEC code 
changes the decay of the relevant two-point correlation 
function from - l/i 2 ( s+1 ' to - 1/^0+1+2^ where n ig 
the number of logical NOTs in a cycle. Such a modifi- 
cation may be overkill for the ohmic case. However, for 
sub-ohmic baths, s < 1, this change in the QEC code 
tremendously improves its effectiveness against correla- 
tions. Of particular importance is the fact that 1// noise, 
so troublesome to qubits based on superconductors, can 
be thought of as the s — > limit of the sub-ohmic case . 
Hence, we have shown that QEC codes can be improved 
so as to be effective in this case. 

The effects of spatial correlations are also straightfor- 
ward to discuss. The form of Eqs. I)ll [l -I)13 [l remains 
the same; the order A 2 terms are simply multiplied by 
an overall prefactor of order 1 which depends on the ge- 
ometry. Thus, applying the logical NOT will reduce the 
effect of spatial correlations between logical qubits sepa- 
rated by a distance larger than v-^A. 



Conclusions — We have shown (1) that the dynamics 
of QEC introduces a natural scale into the problem of 
qubits coupled to a correlated environment, and (2) how 
to use this scale to decrease the computational error. For 
the spin-boson example, the result is a very useful combi- 
nation of QEC with dynamical decoupling: QEC handles 
the short time decoherence while dynamical decoupling 
on the time scale established by QEC reduces the effect 
of long time correlations. 
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